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ABSTRACT 


In  this  report,  a  variational  principle  for  unsteady  body 
wave  problems  is  treated  both  with  and  without  a  convolution 
integral,  and  with  both  linear  and  nonlinear  free  surface  con¬ 
ditions.  Functionals  are  obtained  for  the  numerical  computation 
of  unsteady  flow  fields  near  a  body  that  moves  on  or  beneath  the 
free  surface.  This  formulation  can  be  applied  to  ship  hydro- 
dynamic  performance  problems  of  water  entry  and  body  slamming, 
as  well  as  to  arbitrary  body  motion. 


ADMINISTRATIVE  INFORMATION 

The  work  reported  herein  has  been  supported  by  the  Numerical  Naval  Hydrodynamics 
Program  at  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center.  This 
program  is  jointly  sponsored  by  the  Office  of  Naval  Research  and  DTNSRDC  under  Task 
Area  RR0140302,  Work  Unit  1542-018. 


INTRODUCTION 

In  the  early  1970's  the  David  W.  Taylor  Naval  Ship  Research  and  Development 
Center  (DTNSRDC)  recognized  the  demand  for  advanced  numerical  methods  to  predict  the 
hydrodynamic  performance  characteristics  of  naval  ships,  particularly  when  classical 
methods  proved  inadequate. 

Thus,  in  1974  the  Numerical  Naval  Ship  Hydrodynamics  Program  was  begun  at 
DTNSRDC.  Under  this  program  the  author  previously  investigated  the  steady  ship-wave 
problem  using  a  variational  principle  associated  with  a  localized  finite-element 
technique.  This  method  is  useful  particularly  to  analyze  the  flow  field  near  the 
ship  in  detail;  in  the  far  field,  the  Michell  approximation  can  be  used.  This  report 
extends  the  problem  to  the  unsteady  case. 

For  both  the  steady  and  unsteady  problems,  the  simple  calculation  is  for  a 
linear  free  surface  condition  with  exact  body  boundary  conditions.  An  iterative 
method  is  needed  for  a  nonlinear  free  surface  condition.  However,  in  the  unsteady 
problem,  the  variational  principle  requires  an  integration  with  respect  to  time  using 
the  initial  conditions;  in  this  instance,  a  convolution  integral  is  useful.  The 
variational  principle  for  the  unsteady  body  wave  problem  with  exact  body  boundary 
conditions  is  treated  both  with  and  without  convolution,  and  with  both  linear  and 


*A  complete  listing  of  references  is  given  on  page  13. 


nonlinear  free  surface  conditions.  For  the  linear  free  surface  problem,  a  func¬ 
tional  for  the  variational  principle  is  obtained  with  a  convolution  rather  than  a 
general  integral.  The  convolution  cannot  be  applied  to  a  nonlinear  free  surface 
problem;  the  condition  is  required  for  large  values  of  time.  A  nonlinear  solution 
derived  using  an  iteration  scheme  having  the  linear  convolution  form  is  also  dis¬ 
cussed.  The  time  integration  can  be  eliminated  if  the  motion  is  sinusoidal. 

This  formulation  can  be  applied  to  problems  of  water  entry  and  body  slamming, 
as  well  as  to  arbitrary  body  motion. 

NONLINEAR  PROBLEM 

Since  problems  dealt  with  here  can  be  generalized  easily  to  three-dimensions, 
t or  simplicity  we  first  consider  a  two-dimensional  problem  in  the  rectangular  Car¬ 
tesian  (x,z)  coordinate*  plane.  When  a  body  whose  surface  is  represented  by 


Sg  [ z=h(x, t) ] 


enters  the  water  surface  S„  (z  =  0,  t  <  0;  z  =  h(x,t),  t  >  0)  at  time  t  =  0,  or  when 

r  — 

a  semisubmerged  or  fully  submerged  body  starts  to  move  at  time  t  =  0  and  either  exits 

2 

the  water  or  stops  moving  at  t  =  tj,  then  the  boundary  conditions  for  a  velocity 
potential  4>  are  as  follows: 


4>  =  <l>t  =  0  for  t  <_  0  everywhere 


■|  W)2  -  4>t  +  gh  =  0  on  Sp 


h  -  V<(>V  (h-z)  =  0 


h  -  4)  Kh  +1  =  0 
t  n  x 


on  S  and  S  (t) 
F  s 


Here,  Sg(t)  is  the  submerged  body  surface  varying  with  time  t,  and  n  is  the  normal 
direction  into  the  fluid.  We  consider  potentials  in  the  domain  and  $2  in  t^e 


♦Definition  of  notations  are  given  on  page  iv. 


Then 


domain  D,,,  where  is  the  near  field  including  Ss>  and  is  outside  of  D  . 
at  the  interface  of  and  D^,  we  need  to  have 


d)  =6 
*1  ^2 


In 


-4> 


2n 


(4) 


The  outer  potential  i-n  is  assumed  to  satisfy  the  linear  free  surface  condition 

*2tt  +  s«2z  '  0  <5> 

3  4 

For  such  ^2  we  know  the  time-dependent  Green  function.  ’ 

Now  we  will  construct  a  Lagrangian  for  the  previously  described  problem,  con¬ 
sidering  the  Lagrangian  that  Luke5  used 


J 


dxdt 


where  >>  ,  and  h  vary  with  time,  and  t  is  a  sufficiently  large  time  after  the  body 
has  either  exited  from  or  come  to  rest  in  the  water  so  that  we  can  safelv  assume  that 
the  variation  ant^  |S'^2t  van^3^es  at  t  =  t .  It  will  be  shown  later  that  the 

use  of  a  convolution  integral  necessitates  only  the  initial  condition  without  the 
condition  at  t  =  i . 

Since 


3t 


h(x,t) 

<£  ds  +  <p(z= h) 


(7) 


we  have 
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.T  „  -h(x,t) 


r  h(x.t) 


0  L 


dxdt 


■i  r- : 


X  „T 


dt  <p(z=h )  h 


0  0 


dx 


or 


f  Jj  *t dzd*dt  *  -j  J 

*'«  _  /  \  r\  r» 


0  D1(t) 


h_  4>dxdt  +  J  J  4>dxdz 

0  S  US  D  (t=x) 

s  F  J- 


II 


(8) 


Now  we  Cake  a  variation  of  J  in  Equation  (6),  and  use  the  Green  theorem,  and  the 
condition  at  large  t  =  T  obtaining 


5J 


*J  j  (i  vh  vV’fit+8h 

°  SJ  F 


Sh  dxdt 


-  j"  ("J  V24>1  54)  1  dzdxdt  -  J  J  A’ln’^V*"^  _ 

o  d,  o  slpuss  ' 


4>,  /h“+l  -h  ]  64),  dxdt 
In  x  t  1 


■fi 

0  S, 


-fJ 


(4>in-H>2n)  dzdt  "J  J  (W  6<J>2n  dzdt 

0  S. 


\\  J  <64>2  *2n"*2  **2n>  dzdt  =  ° 


0  S. 


(9) 


However  troni  the  identity. 


4 


{V2J26<t>2-4>2V2(6<l!2)}  dxdz 


=1 


{*2nS<P2 


V*2n>? 


ds  =  0 


'JUS2F 


and  from  the  linear  free  surface  condition  on  Equation  (5),  and  the  condition 

at  t  =  x ,  we  have 


(^V2n'W  dzdt 


(6V2n-W 


dxdt 


=  i J  J  dtdx  =  |  j  (^2i2t^2£<,2t)t-  dx  =  ° 

^2F  ^ 

so  that  the  last  integral  of  Equation  (9)  vanishes.  Since  5$^,  3i> 2 ,  50?n,  and  5h 
are  arbitrary,  we  obtain  from  Equations  (9)  and  (10)  the  corresponding  time- 
dependent,  free  surface  boundary  value  problem  represented  by  Equations  (2)- (4) 
together  with  the  Laplace  equation  for  4> .  Therefore,  solving  the  variational  problem 
with  the  functional  J  of  Equation  (6)  is  equivalent  to  solving  the  Laplace  equation 
with  the  boundary  conditions  as  set  forth  in  Equations  (2)- (4),  provided  we  assume 
that 


6^  =  6<(>2  =  =  0  at  t  =  t 

LINEAR  PROBLEM 

If  we  assume  a  linear  free  surface  condition  in  both  S,_  and  S„_  and  keep  the 

IF  2F 

exact  boundary  condition  [Equation  (3)]  on  the  body  surface. 


>„  •  V 


& 


+1 


(ID 


then  we  may  use 
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J  jj  2  v<|,iV4>i  dzdxdt  ■  Yi  J  {  ^ 


ltt 


dxdt 


0  D 


1 


0  S 


Fo 


/  J  V,  d*dt  -  j  j  (v  I  h)  f2„ 

OS  OS'  / 


dsdt 


0  S  . 


where  S  is  the  projection  of  S.  on  z  =  0.  When  we  take  the  variation  of 
ro  r 


0  J 


"J  II  V\6*l  dzdxdt  _I  I  6^1  (^ln4  i 

0  D,  0  S,„  ' 


*!tt  dxdt 


lFo 


54',  <h  Jh2+l  -h  }  dxdt 
l  in  x  t  / 


0  S 


n 


0  s 


is 


V  &*2n  dzdt  “  I  I  ^ln+V^l  dzdt 


0  S , 


Here,  in  addition  to  Equation  (10),  the  following  equation  holds 


5 1  I  Vltt  dxdt  =J  I  6Vltt  dxdt  +J  J 


itt  dxdt  +  I  !  *i**itt  dxdt 


0  s 


Fo 


0  S 


Fo 


0  S 


Fo 


'  2|  J  «i*m  dxdt  +  j  ^i5^it_^it<S^t)t=T  dx 


0  S 


=  2  J  J  54^4 


ltt 


Fo 


Fo 


0  S 


Fo 


•air  •  . 


(1?) 

J,  we  have 


(12a) 


dxdt 

(13) 
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Thus,  as  in  the  previous  section,  we  can  easily  derive  the  corresponding  linear  free 
surface  boundary  value  problem  using  Equation  (5)  with  the  exact  body  boundary  con¬ 
dition  in  Equation  (11). 


USE  OF  CONVOLUTION 

Equations  (6)  and  (12)  are  Lagrangians  in  a  time-dependent,  two-dimensional 

space  with  nonlinear  and  linear  free  surface  conditions,  respectively.  They  could 

be  localized  in  space  but  not  in  time.  Namely,  we  had  to  specify  the  conditions  on 

S  at  t  =  0  and  t  =  :  with  a  sufficiently  large  t.  On  the  other  hand,  we  did  not 
r 

require  the  initial  condition  *  =  0  at  t  =  0.  In  addition,  such  time  T  when  1 .  =  0 

3  C  1 

on  S  may  be  too  large  for  practical  use.  For  linear  free  surface  boundary  condi- 

tions,  we  can  treat  our  variational  problem  in  the  same  way  as  those  who  have  treated 

variation  principles  for  linear  initial  value  problems^’ ^  using  convolutions  defined 

by 


*  *> 


i 


)1(x,z,t)  4>2(x,z;T-t)  dt 


(14) 


*  v$2 


30,  30-  30, 

nr  *  3^r  +  37“ 


* 


30 


2 


3z 


We  change  Equation  (12)  to 


J 


*  70 ,  dzdx  - 

1  2g 


I  h 


ltt 


dx 


Fo 


0j  dx 


(15) 


If  we  use  the  identity  relation 


and  the  initial  condition  <p^  =  $  ■  0  at  t  *  0,  instead  of  the  conditions  =  0  on 

Sp  at  t  =  0  and  t  =  t  with  large  t  used  in  the  previous  section,  we  obtain,  for  any 
time,  t 


*1  *  *ltt  dX  =j  j  5*l(x* 


t>  ♦ltt^X,T-r^dxdt  + 


J  J  *,<«. 


t)  6<J> ^  (x.T-t)dxdt 


=  j  5fp(x,  t)  f)ltt(x,T-t)dxdt 

0 

Fo 

+  J*  { (x,  t )  5<tlt(x,t-T)-<}>lt(x,t)  6(J>1(x,t-t)}^=0  dx 


=  2  J"  j"  6^-1(x, 


t)  ^ltt (x.T-t)dxdt 


*  htt  dX 


where  i  need  not  be  large,  and  x  represents  a  point  on  the  free  surface  z  =  0.  Thus 
we  obtain  as  in  Equation  (12a) 

6J  =  -  jj  V24>1  *  64>1  dzdx  -  j  6^  *  (<f1ln+  “  <f>lt t  )  dx 
D1  SlFo 

"  J  6*1  *  (♦in  ~  \)dx 


"f 


2)  *  6<D2ndz  -  (1>ln^2n)  *  «Vz 


To  obtain  Equation  (17),  we  used  a  convolution  expression  of  Equation  (16),  where  we 

can  use  the  initial  condition  =  $  =  0. 

Since  6$^,  (S^,  and  5$  are  arbitrary,  we  obtain  the  corresponding  equations 

for  a  time-dependent  linear  free  surface  boundary  value  problem  that  has  a  unique 
3 

solution. 

If  we  lift  out  Sj  so  that  occupies  the  entire  fluid  domain,  then  the  last 
integral  of  Equation  (15)  disappears.  The  resulting  equation  appears  much  simpler 
than  that  obtained  by  Murray2  due  to  a  simple  difference  in  the  treatment  of  the 
free  surface  condition.  Equation  (15)  does  not  give  the  wave  height  as  a  natural 
boundary  condition,  whereas  Murray's  corresponding  equation  does.  However,  from  h  = 
<+> ^ / g  the  wave  height  can  be  obtained. 

If  we  consider  eigen  solutions  that  satisfy  only  the  Laplace  equation;  the 
linear  free  surface  condition.  Equation  (5);  and  the  radiation  condition  in  U  D9, 
then  0,  which  satisfies  the  body  boundary  condition  in  Equation  (11),  can  be  derived 
from  Equations  (15)  and  (17)  by  using 


J 


/h2+l  -  h  ) 
x  t/ 


4>dx 


(18) 


When  we  know  a  functional  whose  minimum  value  is  attained  by  the  solution,  we 

g 

can  find  the  solution  numerically  by  such  methods  as  the  finite-element  technique^ 

9 

or  singularity  method. 

For  example 


<t>  = 


(19) 


where  is  the  Green  function,  which  is  available  for  this  problem  for  a  source 
distribution  on  the  body  surface.  The  source  distribution  nu  will  be  obtained  from 
soltuion  of  the  simultaneous  equations, 


,  1=1,2, 


.N 


(20) 
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SINUSOIDAL  MOTION 


If  we  consider  a  sinusoidal  ship  oscillation  such  as  h  =  f  e^a>t  for  S  of 

t  s 

Equation  (1),  we  substitute 


<f>l  =  ^1(x,z)  e 


icot 


(21) 


into  Equation  (12),  integrate  with  respect  to  t,  and  obtain 


-iwt  , 
e  J 


'  J1  ■  jj  I  v*i»i  -  Tt  J  *1 


dx 


'J  f*l  dx  "J  (Vf*2)^2n  dz 


(22) 


This  is  exactly  the  same  Lagrangian  that  Bai  and  Yeung  used.  Working  from  Equation 
(18),  we  can  apply  eigen  solutions  to  the  whole  field  by  using 


■J  (K-f) 


J^ds 


where 


h  !  =  F  eIu>t 

t  x 

9 

A  similar  functional  was  used  by  Sao  et  al.  to  solve  the  problem  of  a  heaving 
oscillation  of  a  dock. 

ITERATIVE  SCHEME 

For  problems  with  the  linear  free  surface  condition,  we  can  completely  localize 

the  numerical  scheme  in  with  0  ^  t  <  t  for  any  t  with  the  help  of  the  convolution 

form.  However,  for  nonlinear  problems,  the  finite-element  technique  has  to  rely  on 

2 

an  iteration.  We  may  thus  use  an  iterative  free  surface  condition  on  z  =  0 


10 


>itt  +  8  *iz  +  fi  =  0 


ghi  -  *it  +  gi  =  0 


where,  for  the  first  order  perturbation  solution,  f^  =  0,  =  0  and,  for  the  nth 

order,  f  and  are  known  functions  of  of  the  (n-l)th  or  the  lower  order 
solutions. ^ 

Then  the  Lagrangian  for  each  f .  is 


'i  -/jVhi  *v*u 


'll  *  »lltt  dx 


(z=0  )f>D, 


(z=0)HD1 


'H*fl  d‘  +  .fh«**li 


/  (*u"  I  *21)  *  *; 


SJ.  =  0  <25 

1 

where  the  solutions  for  1=1,  2,...n  -  1  should  be  used  to  determine  the  solution 

when  i  =  n.  Equation  (24)  gives  the  wave  height  h  for  each  i. 

If  we  specify  a  time-dependent,  free  surface  pressure  distribution  p  on  the 

proiection  S  of  S  to  z  =  0  instead  of  h  in  S  ,  we  l.jy  use 
K  J  so  s  t  s 

fl  =  =  P/P  on  Sso 


in  Equations  (23)  and  (24)  for  the  first  order,  where  p  is  the  water  density. 


11 


Although  we  have  discussed  the  two-dimensional  time-dependent  problem.  Equations 
(19)  through  (24)  can  be  extended  to  the  three-dimensional  time-dependent  problem. 
That  is 


J 


Ill 


5'*,  * 


V<}>  dxdydz 


1 

2g 


11 

(z=o)nD1 


♦i 


*  ^ltt  dxdy 


-11 


dxdy  + 


(z=o)no 


If. 


■f>  j  dxdy 


{j  (*l~  2*2) 


*  ds 


and  so  on. 


CONCLUDING  REMARKS 

We  have  formed  functionals  with  both  linear  and  nonlinear  free  surface  boundary 
conditions.  For  the  former  but  not  the  latter  case,  we  could  apply  a  convolution 
integral.  However,  the  body  boundary  condition  is  satisfied  exactly  in  both  cases. 

In  many  cases,  the  flow  field  near  an  arbitrary  body  is  of  interest,  and  eigen  solu¬ 
tions  with  linear  free  surface  conditions  are  known.  Even  in  any  large  unsteady 
motion  such  as  ship  slamming,  the  free  surface  condition  for  a  short  period  in  the 
beginning  may  be  linear,  then  the  convolution  may  be  applied  in  the  early  time 
period.  Especially  in  the  slamming  problem,  the  peak  pressure  is  known  to  be  reached 
early  in  the  beginning  and  estimation  of  the  early  pressure  distribution  on  the 
slamming  body  is  required.  With  this  functional,  we  can  find  the  solution  for  an 
arbitrary  body  numerically  by  such  methods  as  the  finite  element  technique  or  singu¬ 
larity  method.  Thus,  a  wide  application  of  such  functionals  can  be  expected. 


12 


REFERENCES 


1.  Yim,  B.,  "A  Variational  Principle  Associated  with  a  Localized  Finite- 
Element  Technique  for  Steady  Ship-Wave  and  Cavity  Problems,"  Proceedings  of  the  First 
International  Conference  on  Numerical  Ship  Hydrodynamics,  David  Taylor  Naval  Ship 
R&D  Center  (1976). 

2.  Wehausen,  J.V.  and  E.V.  Laitone,  "Surface  Waves,"  Encyclopedia  of  Physics, 
Spr inger-Verlag,  Berlin,  Vol.  IX,  pp.  446-778  (1960). 

3.  Stoker,  J.J.,  "Water  Waves,"  Interscience  Publishers,  Inc.,  New  York  (1957), 
pp.  187-196. 

4.  Yim,  B.,  "Investigation  of  Gravity  and  Ventilation  Effects  in  Water  Entry 
of  Thin  Foils,"  Proceedings  of  the  International  Union  of  Theoretical  and  Applied 
Mechanics  Symposium  held  in  Leningrad  (1971),  NAUK  Publishing  House,  Moscow,  pp. 
471-475  (1973). 

5.  Luke,  J.C.,  "A  Variational  Principle  for  a  Fluid  with  a  Free  Surface," 
Journal  of  Fluid  Mechanics,  Vol.  27,  Part  2,  pp.  395-397  (1967). 

6.  Gurtin,  M.E.,  "Variational  Principles  for  Linear  Initial-Value  Problems," 
Quarterly  Journal  of  Applied  Mathematics,  Vol.  22,  pp.  252-256  (1964). 

7.  Murray,  J.C.,  "A  Note  on  Some  Variational  Principles  for  a  Class  of  Linear 
Initial-Boundary  Value  Problems,"  Journal  of  the  Institute  of  Mathematics  and 
Applications,  Vol.  12,  pp.  119-123  (1973). 

8.  Bai,  K.J.  and  R.  Yeung,  "Numerical  Solutions  of  Free-Surface  Flow  Problems," 
10th  Symposium  on  Naval  Hydrodynamics,  Office  of  Naval  Research  (1974). 

9.  Sao,  K.  et  al.,  "On  the  Heaving  Oscillation  of  a  Circular  Dock,"  Journal  of 
the  Society  of  Naval  Architects  of  Japan,  Vol.  130,  pp.  121-131  (1971). 


13 


INITIAL  DISTRIBUTION 


Copies 

1 

2 

4 

2 

1 

2 

1 

1 

13 


12 

I 

I 

1 


CHONR/438  Lee 
NRL 

1  Code  2027 
1  Code  2627 

USNA 

1  Tech  Lib 
1  Nav  Sys  Eng  Dept 
I  Bhattacheryya 
1  Calisal 

NAVPGSCOL 
1  Library 
1  Garrison 

NELC/Lib 

NOSC 

1  Library 
1  Higdon 

NADC 

NCEL/ Code  131 


NAVSEA 


1 

SEA 

031, 

R .  Johnson 

1 

SEA 

031, 

G.  Kerr 

1 

SEA 

031, 

C.  Kennel 

1 

SEA 

05R, 

L.  Benen 

1 

SEA 

05R, 

Dilts 

1 

SEA 

05R, 

N.  Kobitz 

1 

SEA 

05R, 

J.  Schuler 

1 

SEA 

312, 

P . A .  Gale 

1 

SEA 

312, 

J.W.  Kehoe 

l 

SEA 

55W, 

E.N.  Comstock 

1 

SEA 

55W, 

R.G.  Keane,  Jr 

i 

SEA 

61433,  F.  Prout 

1 

PMS 

383, 

Chatterton 

DTIC 

NSF/Engxneering  Lib 


DOT/Lib,  TAD-491.1 
NBS/Rlebanof f 


Copies 

1  MARAD/Lib 

4  U.  of  Cal/Dept  Naval 
Arch,  Berkeley 
1  Eng  Library 
1  Webster 
1  Paulling 
1  Wehausen 

2  U.  of  Cal,  San  Diego 

1  A.T.  Ellis 
1  Scripps  Inst  Lib 

2  CIT 

1  Aero  Lib 
1  T.Y.  Wu 

1  Catholic  U.  of  Amer/ 

Civil  &  Mech  Eng 

1  Colorado  State  U./Eng  Res  Cen 

1  Florida  Atlantic  U. 

1  Tech  Lib 

1  U.  of  Hawaii/St.  Denis 

1  U.  of  Illinois/J.  Robertson 

2  U.  of  Iowa 

1  Library 
1  Landweber 

1  U.  of  Kansas/Civil  Eng  Lib 

1  Lehigh  U. /Fritz  Eng  Lab  Lib 

6  MIT 

1  Yeung 
1  Mandel 
1  Abkowitz 
1  Newman 
1  Sclavounos 
1  Ogilvie 

1  U.  of  Mich/NAME 

1  Library 


15 


PrtECEDLNO  fiiCS  BLANK -NOT  FlUmD 


Copies 

1  U.  of  Notre  Dame 

1  Eng  Lib 

2  New  York  U./Courant  Inst 

1  A.  Peters 
1  J.  Stoker 

4  SIT 

I  Breslin 
1  Savitsky 
1  Dalzell 
l  Kim 

1  U.  of  Texas/Arl  Lib 

2  Southwest  Res  Inst 

1  Applied  Mech  Rev 
1  Abramson 

1  Stanford  Res  Inst/Lib 

2  U.  of  Washington 

1  Eng  Lib 
1  Mech  Eng/Adee 

3  Webb  Inst 

1  Library 
1  Lewis 
1  Ward 

1  Woods  Hole/Ocean  Eng 

1  SNAME/Tech  Lib 

1  Bethlehem  Steel/Sparrows  Point 

1  Bethlehem  Steel/New  York/Lib 

1  Exxon,  NY/Design  Div,  Tank  Dept 

1  General  Dynamics,  EB/Boatwr ight 

1  Gibbs  &  Cox/Tech  Info 

1  Hydronautics/Lib 

1  Lockheed,  Sunnyvale 

1  Potash 


Copies 

1  Newport  News  Shipbuilding/Lib 

1  Ocean ics 

1  Sperry  Rand/Tech  Lib 

1  Sun  Shipbuilding/Chief  Naval  Arch 

2  American  Bureau  of  Shipping 

1  Lib 
1  Cheng 

1  Maritime  Research  Information 

Service 


CENTER  DISTRIBUTION 


Copies 

Code 

Name 

1 

012 

D.  Jewell 

1 

1170 

R.M. 

Stevens 

1 

1170 

G.R. 

Lamb 

1 

1170 

S.  Hawkins 

1 

1500 

W.B. 

Morgan 

1 

1504 

V.  J. 

Monacella 

1 

1520 

w.c. 

Lin 

1 

1521 

w.n. 

Day 

1 

1521 

A.M. 

Reed 

1 

1522 

G.F. 

Dobay 

1 

1522 

M.B. 

Wilson 

1 

1540 

J.H. 

McCarthy 

15 

1540 

B.  Yim 

1 

1542 

T.T. 

Huang 

1 

1542 

J .  Bai 

1 

1560 

D.  Cieslowski 

1 

1560 

Division  Head 

1 

1561 

G.C. 

Cox 

1 

1561 

S.L. 

Bales 

1 

1562 

D.D. 

Mora  ,i 

1 

1562 

E.E. 

Zarnick 

1 

1562 

Y.S. 

Hong 

16 


Code 

1563 

1564 
1564 

5211.1 

522.1 

522.2 


Name 

W.E.  Smith 

J.P.  Feldman 
R.M.  Curphey 

Reports  Distribution 
Library  (C) 

Library  (A) 


17 


I 


